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EXECUTIVE  SUMMARY 


This  report  Illustrates  the  application  of  three  different  procedures  for 
design  of  a  servo-mechanism  controller  which  will  track  input  servo-commands 
and,  at  the  same  time,  accommodate  an  external  disturbance  input.  Two  of  the 
procedures  were  proposed  by  Johnson  [1,  3,  6];  one  uses  optimal  control  tech¬ 
niques  and  the  other  uses  a  linear  algebraic  approach.  The  third  technique, 
for  design  of  a  robust  controller,  was  proposed  by  Davison  [8,  9]. 

These  three  techniques  were  applied  to  an  example  consisting  of  a  linear, 
time-invariant  second-order  plant  with  servo-command  inputs  of  the  form  yc  ■ 

Cq  +  cit  and  an  external  disturbance  input  of  w  ■  efc.  A  servo  controller  was 
designed  using  each  technique  to  demonstrate  the  various  steps  Involved  in 
each  approach.  The  performance  of  the  resulting  controllers  was  then  inves¬ 
tigated. 

The  simplest  controller  design  was  found  to  be  that  resulting  from  the 
application  of  Johnson's  linear  algebraic  approach.  The  controller  designed 
using  Johnson's  optimal  control  approach  was  the  most  computationally  complex 
of  the  three  techniques  and  was  also  the  most  sensitive  to  variations  in  the 
system  gains.  The  controller  resulting  from  Davison's  approach  required  the 
most  additional  Integrators  for  its  Implementation  and  gave  the  worst  tran¬ 
sient  response  to  an  initial  condition.  All  of  the  controllers  were  found  to 
be  sensitive  to  differences  between  the  actual  and  modelled  disturbance  input 
to  the  plant,  with  the  controller  from  Davison's  technique  being  the  least 
sensitive.  All  of  the  controllers  demonstrated  good  servo-tracking  perfor¬ 
mance  in  the  absence  of  external  disturbances  and  all  were  able  to  accommodate 
the  effects  of  the  external  disturbance  out  to  a  run  time  of  about  10  seconds. 
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I.  INTRODUCTION 


Disturbances  are  defined  as  the  uncontrollable  Inputs  which  act  on  a  dy¬ 
namical  system.  There  are  many  varieties  of  disturbance  Inputs  which  can  be 
associated  with  a  controlled  system  and  they  are,  for  the  most  part,  complete¬ 
ly  unpredictable  in  magnitude  and  in  their  arrival  times. 

Johnson  [1-7,10]  introduced  the  idea  of  mathematically  describing  uncer¬ 
tain  waveform-structured  disturbances  by  representing  them  as  a  weighted 
linear  combination  of  known  basis  functions  of  the  form 

n 

w(t)  .VcifiU),  (1) 

i-1 

where  w(t)  is  the  plant  disturbance  vector  and  is  a  p-vector,  and  the 
weighting  coefficients  cj  are  completely  unknown  constants  which  can  change  in 
magnitude  in  a  random,  once-in-a-while  fashion.  The  basis  functions  fj(t)  are 
completely  known  because  they  are  chosen  by  the  designer  based  on  the  waveform 
patterns  exhibited  (or  thought  to  be  exhibited)  by  the  disturbance. 

Johnson  also  proposed  [1,3,6]  two  systemmatic  procedures  for  designing 
multivariable  servomechanism  controllers  which  can  operate  in  the  face  of 
these  unknown,  waveform-type  external  disturbances  and  unknown  waveform-type 
servo-command  inputs.  Davison,  et  al.,  [8,9]  proposed  alternative  design  pro¬ 
cedures  for  the  same  class  of  servomechanism  control  problems. 

In  this  report,  these  design  techniques  will  be  applied  to  an  example 
plant  in  order  to  compare  the  controller  design  procedures.  The  main  steps  in 
each  design  procedure  are  listed.  The  performance  of  each  of  the  controllers 
is  examined  and  the  results  presented.  In  addition,  the  sensitivity  of  each 
of  the  controllers  to  variations  in  the  associated  control  gains  and  to  mis¬ 
matches  between  the  actual  disturbance  acting  on  the  plant  and  the  disturbance 
modelled  in  the  design  process  are  investigated. 


II.  PLANT,  SERVO-COMMAND  AND  EXTERNAL  DISTURBANCE  GENERAL  MODELS 


The  systems  considered  in  this  report  are  described  in  two  different  ways 
to  maintain  consistency  with  the  nomenclature  in  the  references  where  the  de¬ 
sign  techniques  are  presented. 

For  use  in  the  two  design  techniques  proposed  by  Johnson,  these  systems 
are  represented  as  follows.  The  plant  is  described  by  equations  of  the 
general  form 


x(t)  -  Ax(t)  +  Bu(t)  +  Fw(t)  (2) 

y(t)  -  Cx(t)  +  Eu(t)  +  Gw(t)  (3) 

where  x(t)  is  the  plant  state  vector  and  is  an  n-vector,  u(t)  is  the  plant 
control  input  vector  and  is  an  r-vector,  w(t)  is  the  plant  disturbance  vector 
and  is  a  p-vector,  y(t)  is  the  plant  output  vector  and  is  an  m-vector  and  A, 
B,  F,  C,  E  and  G  are  appropriate  size,  known  matrices  with  time-invariant 
elements . 

The  general  form  of  the  disturbance  state  model  is 

w(t)  -  Hz(t)  +  Lx( t)  (4) 

z(t)  -  Dz( t)  +  Mx(t)  +  o(t)  (5) 

where  z(t)  is  the  p-dimenslonal  disturbance  state  vector,  o(t)  is  a  sparsely 
populated  vector  impulse  sequence  and  H,  L,  D  and  M  are  appropriate  size, 
known  matrices. 

The  general  form  of  the  servo-command  state  model  [3]  is 


yc(t)  -  Gcc(t) 

(6) 

c(t)  -  Ecc(t)  +  p(t) 

(7) 

where  c(t)  is  the  v-dlmenslonal  servo-command  state  vector,  p(t)  is  a  sparsely 
populated  vector  Impulse  sequence  and  Gc  and  Ec  are  appropriate  size,  known 
matrices. 

For  use  with  the  design  technique  proposed  by  Davison,  these  systems  are 
represented  as  follows.  The  plant,  in  this  case,  must  be  linear  and  time- 
invariant  and  is  described  as  [8] 

x“Ax  +  Bu  +  E&>  (8) 

y  ■  Cx  +  Du  +  Fu  (9) 

where  u  represents  the  disturbance  vector  and  is  an  O-vector,  x  is  an  n- 
vector,  u  is  an  m-vector  and  y  is  an  r-vector  and  is  the  output  which  is  to  be 
regulated.  The  error  in  the  system  is  the  difference  between  the  output  y  and 
the  specified  reference  input  yref  is  described  as  [9] 


e  »  Cx  +  Du  +  F«i)  -  Gyref . 


(10) 


The  disturbance  vector  is  described  as  [8] 

*1  “  A1  zl» 

«  ■  Ci  z^, 

where  (C^,  A^)  is  observable,  zi(0)  nay  or  may  not  be  known,  and  zj  Is 
vector.  The  specified  reference  Input  vector  is  described  as  [8] 

7ref  "  Ga 

z2  "  ^ 2  z2 
o  -  C2  z2 

where  (C2,  A2)  is  observable,  z2(0)  is  known  and  z2  is  an  n2-vector. 
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III.  MODELS  FOR  THE  EXAMPLE  PLANT 


The  example  [12]  Co  which  each  of  Che  design  cechnlques  Is  Co  be 
applied  is  as  follows*  The  plane  is  described  by 


y  -  y  ■  u  +  w 


The  external  discurbance  which  is  assumed  Co  acC  on  Che  plane  is  described  by 


w  -  w  »  0  , 


almosC  everywhere*  The  servo-command  inpuC  Co  Che  plane  is  given  as 


yc  "  0  . 


almost  everywhere. 


In  terms  of  Che  sCaCe-space  models  given  in  Section  II,  the  differential 
equations  (16),  (17),  and  (18)  can  be  represented  as  follows.  For  Che  plant, 
let 


xi  -  y 


x2  -  y 


Then,  from  Equation  (16),  one  has 


x2»y*y  +  u  +  w"x2+u  +  w 


and  it  follows  chat,  in  terms  of  Equations  (2),  (3),  (8),  and  (9),  the  plant 
model  la  given  by 


/xA  r o  fi  /xA  /o\  +  /o 
l  x2/  1 1  oj  l  x2y  l  i  J  l  i 


y  -  (I,  0)x  . 


The  external  disturbance  given  by  Equation  (17)  represents  an  exponential  type 
disturbance ,  1 *e • , 


w(t)  ■  ce* 


z  *  w  ■  ce‘ 


z-w»,cet-2 


and  Che  external  disturbance  model  can  be  represented  in  terms  of  Equations 
(4)  and  (5)  as 


i  -  z  +  o(t)  .  (2 

In  terms  of  the  disturbance  model  given  by  Equations  (11)  and  (12),  the  dis 
turbance  is  represented  as 


The  servo-command  given  by  Equation  (18)  represents 

y’cCt)  "  C1 

yc(t)  -  cq  +  ext 

which  can  be  expressed  in  terms  of  Equations  (6)  and  (7)  as 

yc  -  (1.  0)c 


G)'  [«  o]  0  +  • 


For  the  model  given  by  Equations  (13)  through  (15),  the  servo-command  is  re 
presented  as 


c  a 


a  -  (1,  0)  /*21\ 


IV.  BACKGROUND 


A  detailed  description  of  the  two  design  techniques  proposed  by  Johnson 
can  be  found  in  [1,3,6].  In  both  procedures,  the  purpose  is  to  design  control 
components,  where  the  total  control  is  allocated  as 

u(t)  -  uc  +  us  ,  (38) 

so  that  the  effects  of  the  external  disturbances  are  absorbed  (or  minimized) 
and  the  primary  control  task  is  achieved.  The  component  u8  is  designed  to 
achieve  the  primary  control  task,  i.e.,  y(t)  ♦  yc(t)  "rapidly.”  The  component 
uc  is  designed  to  accomplish  the  external  disturbance  absorption. 

In  the  procedure  described  in  [1],  optimal  control  techniques  are  applied 
and  the  control  u8  is  designed  to  minimize  a  quadratic  performance  index  J  in¬ 
volving  the  tracking  error,  defined  as 

e(t)  -  yc(t)  -  y(t)  ,  (39) 

and  the  control  component  u8.  In  the  procedure  described  in  [3],  linear  alge¬ 
braic  methods  are  utilized  to  design  u8«  This  technique  is  extended  in  [6] 
where  improved  computational  procedures  are  presented  for  calculating  the 
necessary  feedback  gains. 

Both  of  these  design  techniques  involve  the  design  of  a  composite  state 
observer  to  provide  estimates  of  the  plant,  servo-command  and  external  dis¬ 
turbance  states  for  use  in  the  controllers. 

The  technique  proposed  by  Davison,  et  al.,  [8,9]  is  for  the  design  of  a 
controller  which  is  robust  to  external  disturbance  effects  and  to  perturba¬ 
tions  in  plant  parameters  and  system  gains.  This  controller  requires  the  de¬ 
sign  of  a  stabilizing  compensator,  called  a  complementary  controller,  and  a 
new  type  of  compensator,  called  a  servo-compensator.  The  complementary  con¬ 
troller  is  a  model  of  the  plant  and  is  designed  to  stabilize  the  closed  loop 
system  consisting  of  the  plant/servo-compensator/complementary-controller  com¬ 
bination.  The  servo-compensator  is  a  model  of  the  disturbance/reference  in¬ 
puts  to  the  system.  Its  purpose  is  to  assure  that  the  controlled  system  is 
stabillzable  and  will  achieve  robust  control. 

None  of  the  above  techniques  require  that  the  external  disturbance  be 
measureable.  All  of  the  techniques  use  a  model  of  the  plant  as  part  of  the 
controller.  Davison's  servo-compensator  can  be  designed  without  necessarily 
knowing  the  exact  models  for  the  external  disturbance  or  reference  inputs. 
Johnson's  techniques  also  use  a  model  of  the  external  disturbance  input  and  of 
the  reference  input,  but  these  must  be  either  known  or  assumed  as  known.  In 
Davison's  technique,  consideration  is  restricted  to  linear,  time-invariant 
plants.  Johnson's  techniques  do  not  require  that  the  plant  be  time-invariant. 


V.  DESIGN  TECHNIQUE  1 


la  chis  section,  the  design  technique  proposed  by  Johnson  in  [I]  is  ap¬ 
plied  to  the  plant  of  Section  III.  As  stated  in  Equation  (38),  the  control  u 
is  divided  into  two  parts:  (1)  uc,  which  is  to  be  designed  to  counteract  the 
external  disturbances  and  (2)  us,  which  is  to  be  designed  so  that  y(t)  follows 
yc(t).  To  achieve  the  condition  that  the  output  follows  the  reference  input, 
u8  is  to  be  designed  to  minimize  a  performance  index  J  of  the  form  [1] 

T 

J[u;xo,t0,T]  -  ±eT(T)Se(T)  +  \  J  [  eT(t)Q(t)  e(t)  +  u*  Ru8]dt  ,  (40) 

*=0 

where  S,  Q,  R  are  positive  definite  matrices. 

The  basic  steps  Involved  in  the  design  of  the  control  components  by  this 
technique  are  as  follows: 

1.  Obtain  the  state  model  for  the  expected  external  disturbances  in 
the  form  given  by  Equations  (4)  and  (S). 

2.  Obtain  the  state  model  for  the  expected  servo-commands  in  the  form 
given  by  Equations  (6)  and  (7). 

3.  Check  for  satisfaction  of  the  complete  absorbability  condition  for 
the  external  disturbance. 

4.  If  total  absorption  of  the  external  disturbance  is  possible, 
choose  uc  to  absorb  the  disturbance. 

5.  Implement  a  composite  state  observer  which  will  provide  accurate 
estimates  of  the  plant  and  external  disturbance  states.  Implement  uc  by  using 
the  outputs  of  this  observer. 


6.  Design  u8  by  first  forming  a  composite  "state"  vector 

x  -  (x  |c)T  . 

Next,  express  the  servo-tracking  error  in  terms  of  x  as 

c  -  [-C  |Gc]x  , 

and  rewrite  the  performance  index  given  in  Equation  (40)  as 


1 


«J[«sS3CO*t0»T]  "  2  xt(T)Sx(T)  + 


»/ 


[^(t)Q(t)x(t)-His(t)Ru8(t)]dt 


c0 


(41) 


(42) 


(43) 


7.  Using  standard  linear-quadratic  regulator  theory,  implement  the 
control  u8  as 


ug  "  (Kl<0  |  •  (44) 

Steps  1  and  2  were  accomplished  in  Section  III.  The  disturbance  state 
model  is  given  by  Equations  (27)  and  (28)  and  the  servo-command  state  model  by 
Equations  (33)  and  (34).  For  step  3,  complete  absorption  of  the  external  dis¬ 
turbance  term,  ideally,  will  require  that 


Buc(t)S  -  Fw(t)  (45) 

for  all  admissible  w(t)  and  for  all  t>tQ.  For  this  to  be  possible,  the  fol¬ 
lowing  absorbability  condition  must  be  satisfied: 

R[F]  C  R[B]  ,  to£  t  <  T  ,  (46) 

(where  R[*]  denotes  the  column  range  space  of  [•]),  i.e., 

Rank[B  |  F]  -  Rank[B]  .  (47) 

If  this  condition  is  satisfied,  the  implication  is  that 


F  =  BT 

for  some  matrix  T.  For  this  example. 


Rank[B  F]  -  Rank 


n 


(48) 


(49) 


Rank [ B ]  ■  Rank  [  J -  1  ; 

,1 


(50) 


therefore,  complete  absorption  of  the  external  disturbance  is  possible  and 
P  s  B  T  f  or 


r  -  i  . 

(51) 

The  form  chosen  for  uc  in  [1]  (step  4)  is 

uc  -  Gy  +  HY  , 

(52) 

where  y  is  the  output  of 

V  -  DY  +  Ey  +  (T12  +  DC)Bus 

(53) 

and  is  an  (n+p-m)-vector.  The  terms  shown  in  Equations  (52)  and  (53)  are  de¬ 
fined  in  [1]  as  follows: 


D  -  (T12  +  SC)(AT12  -  T12)  +  T22(DT22  -  i22) 


(54) 


(55) 


E  -  -DE  +  (Tx2  +  BC)(A(C+)t  -  (C+)T)  +  E 


H  *  “IHT22  • 

(56) 

(T  -  -HE 

(57) 

_  T 

T  -1 

T 

T12  “  (T12 

T12  +  T22  T22^ 

T12 

(58) 

_  T 

T  -1 

T 

T22  "  (T12 

T12  +  T22  T22) 

T22  * 

(59) 

The  matrices  T12  sod  T22  **e  a  pair  of  once-dlfferentiable  nx(n+p-m) 
and  mx(tH-p-m)  matrices,  respectively,  such  that 

tc  |  0] [T^2  I  T22]T  =  0  (60) 

and 

T  .  T 

Rank[T12  | T22 ]  =  nfp-rn  .  (61) 

The  matrices  Ti2»  *22.  ^12 »  ^22  are  Part  of  the  structure  of  a  reduced-order 
observer . 

The  reduced-order  observer  (step  5  above)  which  provides  estimates  of 
the  plant  and  disturbance  states  for  use  by  the  controller  Is  given  In  [1]  as 
follows.  Define  the  variable 

e(t)  -  -  5(t)  +  Ejr(t)  +  T12  x(t)  +  T22  z(t)  (62) 

and  let  the  plant  state  estimate  x(t)  and  the  disturbance  state  estimate  z(t) 
be  defined  as 

x(t)  -  [(C+)T  -  T12S}y(t)  +  T125(t)  (63) 

z(t)  -  T22(5(t)  “  Ey(t))  .  (64) 

The  parameter  £(t)  Is  governed  by  the  equation  [1] 

T  •  T 

k  -  (D  +  EHK  +  KT12  +  £C)(A(C+)  -  (C+)  ) 


-  (P  +  EH)E  +  E]y  +  (T12  +  EC)Bu 
-  (D  +  EH) 5  +  1*  +  Qu, 


(65) 


with 


(66) 


a 


a 


«-  (C  | 0] 


A  FH  T12  T12 


t22  t22 


The  evolutloa  equation  for  e(t)  In  Equation  (62)  Is  given  [1]  as 


*e  -  [V  +  EH]  e  +  T22  o(t)  . 


If  £  Is  chosen  In  Equation  (55)  such  that  e(t)  ♦  0  “rapidly,"  then  x  and  z 
will  be  accurate  estimates  of  x  and  w.  This  matrix  £  Is  then  used  In  Equa¬ 
tions  (53),  (54),  (55)  and  (57). 


one  has 


For  this  example,  step  5  continues  as  follows.  From  Equation  (61), 


T  ,  T 


Rank[T12  T^J  -  2+1-1  -  2, 


and  from  Equation  (60)  one  has 


[C  0]  t12 


[1  0 


0]  t12 

t22 


If  Ti2  and  T22  are  chosen  as 


0  0 
0  l 


T22  -  (1.  0)  , 


then  Equation  (70)  Is  satisfied.  From  Equations  (58)  and  (59),  respectively, 
T12  and  T22  can  be  calculated  to  be 


"0  0 
0  1 


t22  -  a,  o)1 


Making  the  appropriate  substitutions  into  Equation  (67)  allows  H  to  be  found 


H  -  (0,  1)  . 


Similarly,  substituting  Into  Equation  (66)  gives  V  as 


1  o" 

e 

_1  0_ 


r't;  •  WV V.  i*.  . 


r;/ .  a* /-V.v.v  v 


*  *  *  >  *  *  *  1  ►  “ 


m 


Using  these  results,  Equation  (68)  can  be  expressed  as 


1  °1  [call  /  T1  oil 

+  (0,1)  >  e  - 

_1  Oj  [012J  )  L1  °12_ 


e  *  Re  • 


It  is  desired  to  have  e(t)  -*•  0  "rapidly,"  therefore,  on  and  012  should  be  de¬ 
signed  to  accomplish  this.  The  characteristic  polynomial  of  R  in  Equation 

(76)  is 

X.2  -  (1  +  oi2)X  -  (on  “  012)  ■  0.  (77) 

For  good  response,  it  is  desired  that  this  characteristic  polynomial  be  given 
by 

X2  +  16X  +  64  -  0  ,  (78) 

l.e.,  that  the  characteristic  roots  be  placed  at  Xi  -  X2  -  -8.  If  Equations 

(77)  and  (78)  are  compared,  it  can  be  seen  that  one  must  have 


•  0  -  0 


The  expression  for  £(t),  Equation  (65),  can  next  be  computed.  The  matrices 
Q  and  T  are  evaluated  to  give 


Q  -  (T12  +  2C)B 


•0- 


?-  (T12  +  EC)[A(C+)t  -  (C+)T]  -  (V  +  EH)  z  +  i  - 


-1296 


where  CT  is  given  by 


C+  -  (CCT)-1C  -(1,0)  , 


and  £  is  finally  expressed  as 


1  -81 

l  - 

1  -17 


/-1296\  (0 

-J  y+  » 


The  plant  and  disturbance  state  estimates  x  and  z  are  then  given  as  shown  in 
Equations  (63)  and  (64),  respectively,  as 

*  -OO 


.o.iv'-.v  ,.-yyy  ■  *y: >  v-  - ’ 
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z 


Cl  +  81y 


(85) 


To  complete  step  5  above,  the  linear  dynamic  device  represented  by 
Equation  (53)  must  also  be  implemented.  From  Equations  (54)  through  (57),  one 
has 


D  - 


E  - 


H  -  (-1,  0) 
G  -  -81  , 


(86) 

(87) 

(88) 
(89) 


therefore,  Equation  (53)  can  be  expressed  as 


(90) 


The  output  of  Equation  (90)  is  then  used,  along  with  the  plant  output  y,  to 
implement  uc,  as  shown  in  Equation  (52),  as 


uc  *  -81y  -  .  (91) 

The  next  step  in  the  design  procedure  is  step  6  above,  the  design  of 
the  control  u8.  In  Equation  (41),  a  composite  ’’state"  vector  was  defined  as 
x  ■  (x  |  c)T,  i.e.,  a  composite  of  the  plant  and  servo-command  states.  The 
performance  index  J,  which  us  is  designed  to  minimize  was  given  by  Equation 
(43)  in  terms  of  this  composite  state  vector.  In  Equation  (43),  the  matrices 

S  and  Q  are  defined  as 


S  -  [-C  |GC]T  S[-C  |gc]  (92) 

Q  -  [-C  |gc]t  Q[-C  |gc]  .  (93) 

For  the  example  in  this  report,  S  and  Q  will  be  chosen  as 

S  -  Q  -  I  .  (94) 

The  gain  matrices  and  K2  shown  in  Equation  (44)  are  given  by  [1,10] 

Ki  -  -R"1  Bt  K*x  (95) 

K2  -  -R"1  BT  K^c  (96) 


12 


•'V-V-VlV-'/J 


and  che  matrices  are  solutions  of  the  coupled  matrix  differential 

equations 


K*  -  (-A  +  BIT1  Bt  Kx)1  Kx  -  K*A  -  CT  Q  C;  Kx(T)  -  CTSC 

KxC  -  (-A  +  BIT1  Bt  KxC)T  KxC  -  KxC  E  +  Ct  Q  Gc;  Kxc(T)— CtSGc 

Equations  (97)  and  (98)  can  be  solved  once  and  for  all  by  integrating 
ward  time,  starting  at  t  ■  T,  using  the  initial  conditions  indicated, 
ing  Equation  (97)  results  in 


fcxl 

*x2 

-,tx3  (kx3“r)-r“rkx2 

• 

• 

-  (1/r) 

2 

_^x3 

kx4_ 

kx4kx3-rkxl-rkx4 

kx4-rkx2-rkx3 

and  expanding  Equation  (98)  results  in 


^xcl  ^xc2 
•  • 

-  (1/r) 

kxc3(kx3“r)+r 

^xc4  ( 'tx3“r )  “rkxc  1 

_'cxc3  kxc4_ 

kx4kxc3~rkXci 

kx4 ^xc4“rkxc2“rkxc3 

For  backward  time  integration,  let 

t  -  T  -  x;  dt  ■  -dx  . 

Substituting  Equations  (101)  into  Equations  (99)  and  (100)  yields 


“<ftxl<*> 

dkx2( t) 

Itx3<kx3“r)-r-rkx2  kx4(kx3~r)-kxir 

dt 

dt 

-  (-1/r) 

dkx3(  t) 

dkx4(  x) 

2 

_kx4kx3”rkxl'"rkx4  kx4~rk*2“rkx3  _ 

dt 

dt 

and 


-dkxd(T)  dkxc2(  -0 
dt  dT 

^kxc3(  dkXc4( 

dt  dt 


(-1/r) 


kxc3(kx3-r)+r 

kx4kxc3_rkxcl 


kxc4(kx3“r)'"rkxcl 


kx4kxc4“rkxc2~rkxc3 


with  the  initial  conditions  given  by 


Kx(T)  -  CTSC  -  CTC 


-1  0 

.  0  0. 


(97) 

(98) 

in  back- 
Expand- 

(99) 


(100) 

(101) 


(102) 


(103) 


(104) 
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KXC(T)  -  -CTSGc 


-CtGc  - 


0 

0 


(105) 


Equations  (102)  and  (103)  were  programmed  In  a  digital  simulation  and 
runs  were  made,  for  several  different  values  of  r,  to  obtain  steady-state 
gains.  The  results  are  presented  in  Table  1. 


TABLE  1.  Steady-State  Optimal  Gains 


R 

k-rl  _ 

kvrl . . 

Xrr4 

10. 

21.231 

20.488 

20.243 

-1.93 

-2.7715 

-0.95345 

-1.8401 

1. 

3.1075 

2.4192 

2.1974 

-1.5538 

-1.7071 

-0.70711 

-1.0987 

0.1 

0.9745 

0.43166 

0.29382 

-0.88591 

-0.4833 

-0.30151 

-0.26711 

0.01 

0.47245 

0.1105 

0.04701 

-0.46777 

-0.11931 

-0.0995 

-0.04654 

0.001 

0.2556 

0.032639 

0.008079 

-0.25537 

-0.0336 

-0.031607 

-0.008071 

To  implement  us,  it  is  necessary  to  provide  estimates  of  the  servo- 
command  states.  This  is  done  by  constructing  an  observer  for  the  servo- 
command  process.  As  given  in  (1],  the  observer  is  constructed  as 


c  -  (E+NG)c  -  Nyc  .  (106) 

The  error  in  the  servo-command  estimate  is  defined  as 

ec  ■  c  -  c  ,  (107) 

and  the  evolution  equation  for  this  can  be  found  as 


ec  -  c  -  c  -  (E+NG)ec  +  tft)  .  (108) 

The  gain  represented  by  the  matrix  N  is  chosen  to  make  ec  0  "rapidly."  For 
this  example,  one  has 


(E+NG) 


\2  - 


Pn  1_ 

Ln2  °J 

and  the  characteristic  polynomial  for  (E+NG)  is  calculated  as 

nl  X  -  n2  ■  0  . 


(109) 


(110) 
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2  -  -8 


It  Is  desired  chat  Che  roots  of  Equation  (110)  be  placed  at  ^  - 
i.e.,  that  the  characteristic  polynomial  be 

2  +  16  +64-0  .  (Ill) 

A  comparison  of  Equations  (110)  and  (111),  results  in 

n^  -  -16,  n2  ■  “64  .  (112) 

Therefore,  Equation  (106)  can  be  represented  as 


and  implemented  to  provide  estimates  of  the  servo-command  states.  The  control 
u8  as  given  by  Equation  (44),  with  and  K2  from  Equations  (95)  and  (96), 
respectively,  can  now  be  implemented  as 

u°  -  -R_1  BT  Kj  x  -  R'1  Bt  Kxc  c 

-  -(l/r)(kx3Xi  +  1^4x2  +  kxc3  ci  +  kxc4  £2)  •  (114) 

Figure  1  is  a  block  diagram  of  the  overall  system.  Including  the 
aervo-command  observer,  the  plant,  the  composite  plant/disturbance  observer, 
and  the  control  components. 

This  system  was  simulated  on  a  digital  computer  and  several  runs  were 
made.  Table  2  is  a  listing  of  the  simulation.  Figure  2  shows  the  response  of 
the  system  to  an  initial  condition  of  y(0)  -  1,  for  each  of  the  values  of  R 
listed  in  Table  1.  Figures  3  and  4  show  the  output  y(t),  and  the  tracking 
error,  respectively,  in  response  to  a  servo-command  input  of  yc  -  1  +  O.lt, 
with  w  -  0,  for  each  value  of  R.  The  best  performance  is  seen  to  be  for  R  - 
0.001.  Figure  5  shows  the  system  performance  in  response  to  the  initial  con¬ 
dition  and  the  external  disturbance  w  -  efc,  with  yc  -  0.  Again,  the  best  per¬ 
formance  is  with  R  -  0.001. 

Figures  6  and  7  show  the  system  output  and  tracking  error,  respec¬ 
tively,  for  a  case  with  R-0.001,  w  -  ec  and  yc  -  1  +  O.lt,  and  Figure  8  gives 
a  plot  of  w  versus  time.  Figures  9,  10,  and  11  are  system  outputs,  with  R  - 
0.001,  to  servo-command  inputs  of  2t,  -2-4t  and  -10+10t,  respectively.  As 
shown,  the  servo-tracking  error  remains  small. 

To  check  the  sensitivity  of  the  controller  performance  to  variations 
in  the  gains,  a  set  of  runs  were  made  with  53!  and  10Z  gain  variations.  Figure 
12  shows  the  servo-tracking  error  for  a  case  with  R  -  0.001,  with  yc  -  0,  w  - 
0,  and  with  system  gain  variations  of  +53!  and  +103!.  Figure  13  shows  similar 
data  for  a  case  with  yc  -  1  +0.1t.  As  shown,  the  system  response  is  very  sen¬ 
sitive  to  small  gain  variations. 


v 


/A 


TABLE  2.  Simulation  Listing  for  System  of  Figure  1 


COMMON  X<8>,DX<S),KUTTA,BT,NX 
DIMENSION  XDAT ( 20 ) 

NVAR-16 

URITE(20)NVAR 


THIS  PROGRAM  SIMULATES  THE  EXAMPLE  PLANT 
IIXTH  THE  CONTROL  DESIGNED  USING 
JOHNSON'S  METHOD  1 


DO  100  1-1,8 
DX(I)*0. 

X<I>-0. 

CONTINUE 
X(4)  «  1. 

TIME-0. 

NX-8 

DT-0.01 

PRINT*/  ENTER  R,UC0,UC1,YC0,YC1  ' 
ACCEPT*, R, IJCO, MCI , YCO, YC1 
PRINT*,'  ENTER  XKX3,XKX4,XKXC3.XKXC4  ' 
ACCEPT* , XKX3, XXX4, XKXC3, XKXC4 
Sll  -  -81. 


S12  »  -17. 

XP1  -  -1296. 

XP2  -  -207. 

EBAR1  -  -1296. 

EBAR2  -  -288. 

XR  -  l./R 
UC  -  0. 

US  •  0. 

U  •  0. 

M  -  0. 

YC  -  0. 

IPRT-0 

CONTINUE 

IF(TIHE.CE.10. )  CO  TO  9999 
IPRT-IPRT+1 
DO  200  KUTTA-1,4 
H  >  WCO*OP(WC1*TIME) 

YC  -  YCO  ♦  YC1*TINE 
SCERR  «  X(2)  -  YC 
5CTK3  »  -SCERR 
DX(1)  •  ^64.*SCERR 
DX(2)  -  XC1)  -  16.*SCERR 


US  -  -XR*<XKXC4*X<1)*XKXC3*X(2}*XKX3*X<4)*XKX4*(-S12*X(4)*X<S>>) 
UC  •  Sll*X(4)  -  X(7) 

U  •  US  ♦  UC 
0X(3)  »  U  ♦  W  ♦  X<4> 

DX(4)  -  X(3) 

DX(5)  *  XP1*X<4)  ♦  S11*X(6)  ♦  X(S) 

DX(6>  •  X(6)  ♦  U  ♦  XP2*X(4)  ♦  612*X(6> 

X2HAT  •  -S12»X(4>  ♦  X(6) 

0X17)  •  EBAR1»X(4)  ♦  X(7)  ♦  S11*X(8) 

DX(8)  *  ESAR2*X(4)  ♦  US  ♦  S12«X<8) 

GO  TO  <30,60, 30,40).KUTTA 

CONTINUE 
TIME-TIME*. 5*DT 
CONTINUE 
CALL  RUNE 


. V'*’  '• 

a.’ A  *, 


ISwwh 


♦  ^ •  ,  * 


V  -  ’ .  V  %  ‘  „  V* 


TABLE  2.  Simulation  Listing  for  System  of  Figure  1  -  Continued 


200  CONTINUE 

TRKERR  ■  YC  -  X<4> 

EXER  -  W  -  X<7) 

XMTdMU) 

XDAT(2)-X(2> 

XDAT(3)-X(3> 

XDAT(4)«X<4) 

XSAT(B)«X<B) 

XDAT(6)*X(£) 

XBAT(7)-X<7> 

XDAT<8>«XC8> 

XDAT(9)-YC 
XDAT  <  10 )-TIHC 
XMT(il)>U8 
XBAT (12) >UC 
XDAT(13)-X2HAT 
XDAT(14)*TRKERR 
XBAT(19)«SCTKE 
X0AT(i6)m 

IF(IPRT.NE.IO)  CO  TO  BOO 
NRITEC20)  (XDAT(I),I*1,MVAR) 

IPRT-0 

_  «T>  -  ' ,X(4), '  YC  -  ',YC 

900  CO  TO  1000 

99  CONTINUE 
STOP 
END 

SUBROUTINE  RUNK 

COHNON  X(8),DX(8),KUTTA,DT,NX 


10 


20 

30 


40 

90 

60 

70 

80 


DIMENSION  XA(8),DXA(8) 

CO  TO  (10,30,50,70), KUTTA 
DO  20  I-1,NX 
XA(I)-X(I) 

DXA(I>*DT«DX(I) 

X(I)“X(I)».S«DXA(I) 

RETURN 

TBT-2.0DT 

HDT-.5*DT 

DO  40  I«1,NX 

DXA(X)*DXA(I)+T9T«DX(X) 

X(I)-XAU)*HDT»DX<I) 

RETURN 

DO  60  I-l.NX 

VDT-DT»DX(I> 

DXA(I)-DXA(I)+2.*m 

X(I)«XA(I)*VDT 

RETURN 

DO  80  1*1, MX 

X(I)*XA(I)+(DXA(I)*DT*DX(I))/6. 


RETURN 

END 


TRACKING  ERROR 


TIME  <  SEC  > 


TIME  <SEC> 


Figure  12.  Servo-tracking  error,  with  R  m  0.001,  with  system  gain 
variations,  for  a  case  with  yc  ■  0,  w  *  0,  y<0)  *  1. 


Figure  13.  Servo- tracking  error,  with  R  ■  0.001,  with  system  gain 
variations,  for  a  case  with  yc  -  1  +  O.lt,  v  *  eS 
y(Q)  «  1.  _ 
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VI.  DESIGN  TECHNIQUE  2 


In  this  section,  the  design  technique  proposed  by  Johnson  in  [3,6]  will  be 
applied  to  the  plant  of  Section  III.  As  was  the  case  in  the  procedure  of 
Section  V,  the  control  u  is  divided  into  two  parts,  and  u8.  These  two  con¬ 

trol  components  have  the  same  tasks  for  the  design  of  this  section  as  they 
had  in  the  previous  technique;  however,  the  form  chosen  for  uc  is  different 
and  Ug  will  be  designed  by  using  linear  algebraic  techniques. 

The  basic  steps  involved  in  this  technique  are  as  follows: 

1.  Obtain  the  state  model  for  the  expected  external  disturbances  in 
Che  form  given  by  Equations  (4)  and  (5). 

2.  Obtain  the  state  model  for  the  expected  servo-commands  in  the  form 
given  by  Equations  (6)  and  (7). 

3.  Check  to  see  if  the  effect  of  the  disturbances  on  the  output  can 
be  completely  absorbed,  either  by  choosing 

uc  -  Az  (115) 

such  that 

BA  +  FH  =0  ,  (116) 

or 

C(t)«(t,T)[B(T)A(T)  +  F(t)H(x)J  50,  (117) 

i.e.,  (B A  +  FH)  is  in  the  null  space  of  C(t)*(t,T),  in  which  case  the  disturb¬ 
ance  effects  will  be  unobservable  in  the  output  space. 

4.  To  design  u8,  first  check  for  satisfaction  of  the  "exact  track- 
ability"  condition  [3]  for  exact  servo- tracking ,  i.e., 

R[Gc(t)]  9  R[C(t)J  ,  t0  <  t  <  T  .  (118) 

If  Equation  (118)  is  satisfied,  then  it  is  possible  to  express  Gc(t)  as 

Gc(t)  -  C(t)0(t)  ,  t0  <  t  <  T  (119) 

for  some  (possibly  non-unique)  matrix  0(t)«  If  this  condition  is  satisfied, 
then  from  the  expression  for  servo  error  ey,  defined  as 

ey  ■  yc  “  y  ■  Gcc  -  Cx  -  Gc(9c  —  x)  ,  (120) 

a  new  variable  e88  is  defined  as  [3] 

e88  -  0c  -  x  .  (121) 

5.  The  control  us  is  of  the  form  [3] 

us  "  Si(t)x  +  S2(t)c  .  (122) 
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Using  this  definition  for  u8,  the  evolution  equation  for  eas  can  be  derived 
as  [3] 

e88  “  (A+BSi)e89  +  [0E  -  A0  +  9  -  B(Si©  +  S2>]c  -  Bz  +  Q|i(t) 

-  (A  +  BSi)ess  -  Vc  -  Bz  +  0p(t)  ,  (123) 

where  V  is  defined  as 


V  --0E  +  A0-0  +  B(Si©  +  S2)  •  (124) 

Then,  us  (i.e.,  S^,  S2)  is  designed  so  that  ess  rapidly  approaches  the  null 
space  of  C.  This  design  procedure  is  as  follows: 

a.  Si  is  designed  so  that  for  the  homogeneous  differential  equation 
associated  with  the  e88  expression,  e88  +  N(C). 

b.  S2  is  designed  so  that  the  remaining  terms  in  the  es8  expression 
which  can  affect  e8S  are  either  zeroed  out  or  are  confined  to  the  null  space 
of  C*.  The  terms  in  question  are  given  by  Vc,  as  defined  in  Equations  (110) 

and  (111),  since  Bz  is  in  the  null  space  of  C$  and  will  not  affect  the  output. 
Thus,  S2  is  designed  so  that 

C(  t)  ®(t ,  x)  [  0(  x)E(  x)-A(  x)  0(  x)+0(  x)  -  B(T)Si(r)0(T)  -  B(x)S2(x)] 

-  C(t)«(t,x)V(x)  =  0  .  (125) 

6.  Finally,  the  state  reconstructors  are  designed  to  give  estimates 
of  the  plant,  disturbance  and  servo-command  states  to  permit  practical  imple¬ 
mentation  of  uc  and  u8. 

Steps  1  and  2  above,  were  accomplished  in  Section  III,  with  the 
disturbance  state  model  given  by  Equations  (27)  and  (28)  and  the  servo-command 
model  by  Equations  (33)  and  (34).  For  step  3,  since  the  plant  in  this  example 
is  time-invariant,  the  requirement  given  by  Equation  (117)  reduces  to  the  re¬ 
quirement  [3]  that 

C[B  |AB]  -  0,  (126) 

where 

B  -  BA  +  FH  .  (127) 

Making  appropriate  substitutions  into  Equation  (127)  results  in 

•-CK) 
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and  it  follows  that 


A  -  -1  (129) 

will  satisfy  both  Equation  (116),  and  Equation  (126).  The  control  uc,  there¬ 
fore,  is  chosen  as 


,  uc  »  Az  “  -z  .  (130) 

I 

In  step  4,  a  check  must  be  made  for  satisfaction  of  the  exact  track- 
ability  condition.  For  this  example. 


rank[Gc  |C]  -  rank[l  0  |  1  0]  -  1 
rank[Gc]  -  1, 


(131) 


hence.  Equation  (118)  is  satisfied,  and  from  Equation  (119)  it  is  found  that 


0-1 


(132) 


To  proceed  with  step  5,  the  design  of  the  control  u3,  from  step  5. a. 
above,  the  gain  matrix  is  to  be  designed  so  that  the  solution  of  the  homog¬ 
eneous  differential  equation 


ess  -  (A  +  BS!)ess  (133) 

approaches  N(C)  rapidly.  For  this  example,  Equation  (133)  can  be  expanded  as 


-ss 


0  1 
l+su  *12 


ess  ■  Aless 


(134) 


and  the  characteristic  polynomial  of  Ai  is  found  to  be 

X2  -  si2k  -  (sn  +  1)  -  0  .  (135) 

For  good  response,  it  is  desired  that  the  characteristic  roots  of  Aj  be  lo¬ 
cated  at 

Al  -  X2  -  -10  ,  (136) 

i.e.,  that  the  characteristic  polynomial  of  Ai  be  given  by 

X2  +  20  X  +  100  -  0  .  (137) 

^  Equations  (137)  and  (135)  will  be  equal  if 

su  -  -101,  si2  -  -20  .  (138) 

i 
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w  .  v  .  ^  w .  r;  ^  r.  ^  .  r,  *  ,  *  .  .  •  ,  *  .  *  ,  •  .  • 


The  next  step  in  the  design  procedure  is  step  5.b.,  the  design  of  S2. 
As  stated  earlier,  S2  must  be  designed  so  that  Equation  (125)  is  satisfied. 
From  [3],  in  order  to  satisfy  Equation  (125)  it  is  necessary  and  sufficient 
that  an  S2  exist  so  that 

BS2  =  V  +  0E-A0+  0  -  BSiQ,  t0<  t  <  T  (139) 

for  some  matrix  V  which  satisfies  C(t)<I(t,T)  V(  t)  =0,  tQ  t  t,  tg  £  t  T. 

Therefore,  find  a  V  and  then  solve  for  S2  from  Equation  (139).  To  satisfy 
Equation  (139),  V  must  satisfy  the  condition  [3] 

rank[B  |  V  +  0E-A0  +  0]  -  rank[B]  .  (140) 

If  Equation  (140)  is  satisfied,  then  [3] 

V  +  GE  -  A0  +  0  ■  BE  (141) 


for  some,  possibly  non-unique,  matrix  E. 

For  this  example.  Equation  (140)  is  expanded  as 


rank[B  V  +  ®-A0  +  0]  -  rank 


rank[B]  -  rank 

and  it  can  be  seen  that 

rank 


0 


V  + 


0  0 
-1  0 


1  , 


(142) 


(143) 


V  + 


:] 


-  1 


(144) 


if  V  ia  chosen  to  be  0.  This  choice  for  V  will  also  satisfy  Equation  (125). 
Given  this  choice  for  V,  Equation  (141)  can  be  expanded  as 


0  0  0  0 
|-1  0  qi  0^2 

and  therefore  I  is  given  by 

E  -  (-1,  0)  . 

Next,  find  an  S2  which  will  satisfy  Equation  (139)  by  setting  [3] 


(145) 


(146) 


S2  -  E  -  S !0  .  (147) 

From  Equation  (147)  then,  S2  is  found  as 


S2  -  (-1,  0)  -  (-101,  -20)  -  (100,  20)  . 


(148) 


The  control  us,  as  given  by  Equation  (122),  can  now  be  expressed  as 

u8  -  S]X  +  S2c  -  (-101,  -20)x  +  (100,  20)c  ,  (149) 

and  the  total  control  vector  u  Is 

u  •  uc  +  ug  -  -z  -lOlxi  -20x2  +100ci  +20c2  •  (150) 

For  practical  Implementation,  estimates  of  the  plant,  servo-command,  and 
disturbance  states  must  be  used  and  Equation  (150)  Implemented  as 

u  ■  -z  -  lOlxi  -  20x2  +  100ci  +  20c2  .  (151) 

The  estimates  of  the  plant  and  disturbance  states  will  be  provided  by 
a  full-order  composite  state  observer  designed  using  a  "recipe'*  provided  In 
[3].  The  observer  Is  defined  by 


(152) 


By  defining  the  estimation  error  as  (  e*  |  -  (x  |  z)T  -  (x  j  z)T,  the  error 

between  Impulses  of  o(t)  can  be  shown  to  obey  the  homogeneous  differential 
equation 


The  gain  matrices  K^,  K2  are  designed  so  that  (  e*  |  tz)^  >0  "rapidly." 

Using  the  plant/disturbance  models  given  In  Section  III,  Equation 
(153)  Is  expanded  to  give 


The  characteristic  polynomial  of  A2  Is 

X3  -  (kn  +1)X2  +  (kn  -k12  -1)X  +  (l+k12  -k21)  -  0  . 
It  Is  desired  that  the  observer  poles  be  located  at 

Xi  •  -4,  X2  -  -5,  X3  -  -7  , 
which  would  result  In  a  characteristic  polynomial  of 


(154) 

(155) 

(156) 
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(a. 

”-17  1 

-100  0 

o' 

1 
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The  observer  used  for  the  servo-command  state  estimates  is  the  same  one  de¬ 
signed  in  Section  V,  Equation  (113). 

Figure  16  shows  a  block  diagram  of  the  plant/observer/controller  com¬ 
posite  system  for  the  design  of  this  Section.  This  system  was  simulated  on  a 
digital  computer  and  several  runs  were  made  to  check  the  performance  of  this 
system  (see  Table  3  for  a  simulation  listing). 

Figure  17  gives  the  tracking  error  for  two  cases:  (1)  yc  -  0,  w  -  0 
and  (2)  yc  ■  0,  w  -  efc.  The  curve  for  case  (1)  corresponds  to  the  results 
shown  in  Figure  2  for  the  controller  of  Section  V,  and  seems  to  yield  about 
the  same  response  as  the  case  with  R  -  0.001  on  Figure  2.  Figures  18  and  19 
show  the  output  and  the  tracking  error,  respectively,  for  a  case  with  yc  - 
1  +  O.lt,  w  ■  0.  These  results  correspond  to  those  of  Figures  3  and  4.  The 
settling  times  of  the  transients  agree  with  the  results  in  Figures  3  and  4  for 
the  case  with  R  -  0.001,  but  the  magnitudes  of  the  excursions  are  less  for  the 
controller  of  this  Section.  Figure  20  shows  the  output  response  for  a  case 
with  yc  “  1  +  O.lt,  w  ■  ec,  corresponding  to  the  results  shown  in  Figure  6. 

The  transients  in  this  case  are  again  smaller  in  magnitude  than  those  in 
Figure  6.  Figure  21  presents  the  tracking  errors  for  four  cases:  (1)  yc  ■ 

1  +  O.lt;  (2)  yc  -  2t;  (3)  yc  ■  -2-4t;  (4)  yc  -  -10+10t.  These  curves  cor¬ 
respond  to  the  data  shown  in  Figures  7,  9,  10,  and  11,  respectively.  The  re¬ 
sults  are  equivalent. 

Several  runs  were  made  to  check  the  sensitivity  of  the  controller  de¬ 
sign  of  this  section  to  variations  in  the  gains.  Figure  22  is  a  repeat  of  the 
case  in  Figure  16  with  yc  •  0,  w  -  0,  but  with  +52  and  +10Z  variations  in  the 
gains  shown  on  Figure  16.  The  results  correspond  to  the  results  in  Figure  12 
for  the  design  in  Section  V.  As  shown,  the  system  in  this  section  is  not  as 
sensitive,  even  to  +10Z  gain  variations,  as  was  the  system  of  Section  V.  In 
fact,  as  shown  in  Figures  23  and  24,  the  controller  of  this  section  can  with¬ 
stand  gain  variations  from  -25Z  to  +50Z  without  significant  performance  degra¬ 
dation. 


Cent  block  diagram,  technique 


TABLE  3.  Simulation  Listing  for  System  of  Figure  16 


X*7),MU7>,iaJTTA,BT,IIX 
DIMENSION  XDAT(20) 

MVMt-lfi 
URITE(20)NVAR 

THIS  PROCRAH  SSSSTu  -  EXP(TIHE) 
WITH  A  STATE  OBSERVER  FOR  A  SERVO-PROBLEM 
FOR  JOHNSON'S  METHOD  2 


* 

* 


DO  100  1*1,7 

DX(I)*0. 

X(I)*0. 

100  CONTINUE 
X(4)  *  1. 
TIME-0. 

NX-7 


DT-0.01 

PRINT*,'  ENTER  HC0,UC1,YC0,YC1  ' 
ACCEPT*, UC0,HC1,YC0,YC1 
PRINT*, '  ENTER  CAIN  MULTIPLIER  ' 
ACCEPT*,  OPR 
mi  -  -16.*<1.*CHPR> 


XN2  «  -S4.*(1.«CMPR) 

511  -  -10i.*U.*CNPR> 

512  -  -20.*<1.*CHPR> 

XTll  -  -t7.*<l.*CNPRJ 
XK12  «  -101.*(1.«CMPR) 

XK21  -  -240.*<1.*CHPR> 

521  *  100.*<1.*CHPR> 

522  «  20.X1.4CMPR) 

UC  •  0. 

US  *  0. 

U  *  0. 

U  •  0. 

YC  -  0. 

IPRT-0 

1000  CONTINUE 

IFITIHE.CE.IO.)  CO  TO  9999 
IPRT-IPRT+1 
DO  200  KUTTA-1,4 
U  *  HCO*EXP(UCl*riME) 

TC  -  YCO  ♦  YC1»TIHE 
SCERR  -  X(2)  -  YC 
BCTKE  -  -SCERR 


DX(1)  *  -64.*SCERR 
DX(2)  -  XC1)  -  1S.*SCERR 

US  -  S11*X(4>  ♦  S12*X(6)  ♦  S21*X(2)  ♦  S22*X(1> 
UC  •  -X<7) 

U  -  US  *  UC 

DX(3)  »  U  ♦  H  ♦  X(4) 

DX(4)  -  XC3) 

DX<5>  *  -XK11»X(4)  ♦  XK11*X(9)  ♦  X(6) 

DX(S>  •  (l.*XK12)*X(3)  -  XK12*X(4)  ♦  US 
DX(7)  -  -XJC21*X<4)  ♦  XK21*X<8)  ♦  X(7> 

CO  TO  <30,60.30,40>,KUTTA 
30  CONTINUE 

TIME-TIME*. 6*DT 
40  CONTINUE 
SO  CALL  RUNK 
200  CONTINUE 

TRXERR  -  YC  -  X<4> 

EXER  -  U  -  X(7) 

XDAT(1)-X(1) 

XDAT(2)-X(2> 

XDAT(3)-X(3) 

XDAT(4)-X(4> 


uu 


TABLE  3.  Simulation  Listing  for  System  of  Figure  16  -  Continued 


900 
9999 


XB*T<«)aX<C» 

XMT(S)"X<6) 

XBAT(7}*X<7) 

XBAT(8)*£XER 

XBAT<9)*TC 

XDAT(10)*TIHE 

XDAT<ll)*US 

XBAT(12)*UC 

XMT(13)«X2HAT 

XBAT<14)*TRrERR 

XMT(19)aSCTKE 

XMT(16)-U 


IFUPRT.NE.10)  CO  TO  900 
HRITE(20)  <X9AT(I),I*1,NVAR) 

IPRT-0 

pjnjnra,'  tine  ■  mine,'  rm  -  \xu>,'  rc  - 

JC  ■  *.TC.'  TCHAT  ■  ',X(2), '  U  » 
X2HAT  -  '.X2HAT,'  X<2>  - 

CO  TO  1000 

CONTINUE 

STOP 


10 

20 

30 


40 

90 

60 

70 

80 


SUBROUTINE  RUNK 

CONNOR  X(7),DX(7),KUTTA,DT,NX 

DIMENSION  XA(7),DXA(7) 

«*  TO  <10, 30 ,90 ,70), run* 

00  20  1*1, NX 
XAU)-X(I) 

BX*< I )*DT*DX( I ) 
X(I)-X(I)+.9»PXA(I) 

RETURN 

TDT*2.«DT 

HDT-.9#DT 

BO  40  1*1, NX 

DXA<I)-DXA(I)*TDT«DX<I) 

X(I)*XA<I)+HDT»DX(I) 

RETURN 
BO  60  1*1, NX 
VBT*0T«BXU) 

BXA( I )*DXA( I ) *2. «VBT 

X<I)-XA<IHVPT 

RETURN 


00  80  1*1, NX 

X<I)«XA<I)4<0XA(I)4DT*DX(l))/6. 


'.rc 

u 

'»X(3> 


TRACKING  ERROR  «t  TRACKING  ERROR 
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Two  runs  were  made  to  Investigate  the  robustness  of  this  controller  to 
differences  between  the  actual  and  assumed  external  disturbance  inputs  to  the 
plant.  Figures  25  and  26  present  the  results,  versus  a  reference  curve  for  w 
-  ec,  for  cases  with  w  -  e°*5t  and  w  »  e1*5*,  respectively.  As  shown,  this 
controller  is  also  sensitive  to  variations  from  the  assumed  disturbance  model. 


VII.  DESIGN  TECHNIQUE  3 


In  this  section,  the  design  technique  proposed  by  Davison  and  Goldenberg 
In  [8]  and  by  Davison  In  [9]  for  design  of  a  "robust  feedback  controller"  is 
applied  to  the  design  of  a  controller  for  the  plant  of  Section  III.  The  ro¬ 
bust  feedback  controller  consists  of  [8]  a  "servo-compensator"  and  a  "stabil¬ 
izing  compensator."  This  design  technique  assumes  that  the  plant  is  linear 
and  time-invariant  and  is  described  by  the  general  form  given  in  Equations  (8) 
and  (9).  The  tracking  error  is  the  difference  between  the  output  and  the 
specified  reference  input  and  is  given  by  Equation  (10).  The  general  form  of 
the  model  for  the  external  disturbance  is  given  by  Equations  (11)  and  (12)  and 
the  general  form  for  the  servo-command  input  is  given  by  Equations  (13),  (14) 
and  (15) . 

The  steps  which  comprise  the  design  procedure  for  the  robust  feedback 
controller  are  as  follows: 

1.  Check  for  satisfaction  of  the  necessary  and  sufficient  condition 
for  the  existence  of  a  robust  controller  for  the  system.  In  order  that  a 
robust  controller  exist,  the  following  conditions  must  all  hold  [8]: 

a.  (A,B)  must  be  stabllizable. 

b.  (C^A)  must  be  detectable,  where  ym  -  C^x  +  D,,^  +  Fm  u>  are 
the  only  outputs  which  are  available  for  measurement. 


m  >  r 


(160) 


d.  The  transmission  zeros  of  (C,A,B,D)  do  not  coincide  with  Xj, 

i  •  1,2, - ,q • 

e.  ym  must  contain  the  actual  output  y. 

Conditions  c  and  d  are  equivalent  to  satisfaction  of  [8] 


A  -  *4  I  B" 

■  n 

C  D 


+r,i—l,2,  » q  ■ 


(161) 


To  perform  this  step,  it  is  necessary  to  have  the  minimal  polynomials  of  the 
matrices  used  to  model  the  disturbances  and  servo-commands  (A}  and  A2  in 
Equations  (11)  and  (14)).  The  X*  are  the  zeros  of  the  least  common  multiple 
of  these  two  minimal  polynomials. 

The  form  of  the  controller,  if  one  exists,  is  chosen  as  [8] 


u  ■  Kq  x  +  , 


(162) 


where  £  is  the  "servo-compensator"  output  and  is  an  rq-vector  and  x  is  the 
"stabilizing  compensator"  output. 


2.  Design  the  servo-compensator,  according  to  a  recipe  given  in  [8]. 
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3.  Since  Che  plane  in  this  example  is  unstable,  a  full-order  stabil¬ 
izing  observer  is  designed  and  applied  to  it  to  provide  a  stable,  satisfactory 
dynamic  response. 

4.  The  controller  gains  Kq  and  K  in  Equation  (162)  are  chosen  so  that 
the  composite  system  has  a  satisfactory  dynamic  response  and  the  complementary 
controller  is  implemented  to  provide  the  estimates  of  the  plant  states  used  in 
the  control  u. 


To  proceed  with  step  1,  first  determine  the  minimal  polynomials  of 
and  A2.  The  matrix  \  is  given  in  Equation  (29)  as  ^  ■  1<  The  minimal  poly¬ 
nomial  of  A^,  Ax,  is  found  as 


|  XI  -  Ax  |  -  X  -  1  , 

Ax  -  X  -  1  . 

The  matrix  A2  is  given  in  Equation  (35)  as 


(163) 

(164) 


(165) 


and  the  minimal  polynomial  of  A2,  A2,  is  found  as  follows.  The  characteri¬ 
stic  polynomial  of  A2  is 


XI  -  A2 


X  -1 
0  X 


\2  . 


(166) 


Since  A2  is  not  a  null  matrix,  but  ^2  is  a  null  matrix,  the  minimal  polynomial 
of  A2  is  given  by  ^ 

A2  -  X2  .  (167) 

The  least  common  multiple  of  Ax  and  A2  is 

A  -  X2  (X  -  1)  ,  (168) 


and  the  zeros  are  seen  to  be 


Xx  -  0,  X2  -  0,  X3  -  1  .  (169) 

Proceeding  with  step  l.a.,  (A,B)  must  be  checked  for  controlla¬ 
bility.  The  pair  (A,B)  will  be  completely  state  controllable  if,  and  only  if, 
the  composite  matrix  Pi,  where 


Pi  -  [B|AB]  ,  (170) 

is  of  rank  n  [11].  For  this  example  then,  with  A  and  B  as  given  in  Equation 

(22), 
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rank  [PjJ  ■  2  *  n  , 

therefore,  the  pair  (A,B)  is  completely  state  controllable. 


(171) 


(172) 


From  step  l.b.,  the  pair  (C, A)  must  be  checked  for  observability. 
The  pair  (C,A)  is  completely  observable  if,  and  only  if,  the  composite  matrix 
P2,  where 

P2  -  [CT  I  At  Ct1  ,  (173) 


(174) 


P2  -  [C*  |  AT  Ct]  , 

is  of  rank  n  [11].  For  this  example,  with  C  given  in  Equation  (23), 


1  0 

m  9 

0  1 


rank[P2]  ■  2  *  n  , 

therefore,  the  pair  (C,A)  is  completely  observable. 


(175) 


To  check  for  satisfaction  of  the  conditions  in  steps  l.c.  and 
l.d..  Equation  (161)  will  be  used  with  the  X's  of  Equation  (169).  For  X^  *  0, 
the  result  is 


0  1  0 

rank  1 _ 0__1_  «3-n  +  r»3  . 

10  0 


(176) 


For  >2  “  0»  the  result  is  the  same  as  in  Equation  (176).  For  X3  ■  1,  the  re¬ 
sult  is 


-110 
1-11 
1  010' 


n  +  r  *  3  . 


(177) 


Equation  (161)  is  thus  satisfied  for  each  X. 

Since  y,  the  actual  output,  is  assumed  to  be  measureable,  step 
l.e.  is  satisfied.  Therefore,  a  robust  controller  exists  for  this  system. 

The  next  step  in  the  design  procedure  is  2.,  the  design  of  the 
servo-compensator.  The  general  servo-compensator  is  described  by  [8] 


l  -  C*£  +  B*e, 

where  e  is  given  by  Equation  (10),  and  [8] 

B*  -  tB 


(178) 


(179) 
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C*  -  x  block  dlag(C,  C, 


— ,  C)t_1 


(180) 


r  matrices 


with  x  a  non-singular  real  matrix;  8  a  real  matrix  of  rank  r  such  that  (block 

diag  (C,  C, - ,  C),  8 }  is  controllable;  and  C  is  defined  as  a  qxq  companion 

matrix  of  the  form  [8] 


0  10 - 0 

0  0  1 - 0 

II!  II 

L  I  I  II 

III  II 

III  II 

-Si  -&2  -63 - 6, 


(181) 


The  coefficients  5^ ,  63, 
the  polynomial 


■,  in  Equation  (181)  are  the  coefficients  of 


6  xq_1  + - +  62  x  +  61,  -  n  (x  -  Xi), 

q  i-1 


(182) 


where  the  Xj  are  the  zeros  of  A.  Expanding  Equation  (182)  using  the  X's  of 
Equation  (169)  results  in 


(X  -  Xi)(X  -  X2)(X  -  X3)  -  X2  (X  -  1)  -  X3  -  X2  , 
and  the  coefficients  are  seen  to  be 

61  -  0,  62  -  0,  63  -  -1  . 

For  this  example  then,  the  matrix  C  of  Equation  (181)  is  given  by 


0  10 

-001. 

0  0  1 


(183) 


( 184 ) 


(185) 


Following  a  suggestion  in  [8],  let  the  matrix  x  -  I  and  let  B  - 

block  diag  (yi.  Y2, - ,  yr)  where  yi  -  Y2  “ - -  yr  -  (0  0 - 0  1)T. 

In  this  case,  since  r  -  1,  one  has 


8*  -  -cB  -  IB  -  (0  0  1)T  , 


(186) 


0  10 

-001. 

0  0  1 


(187) 


The  servo-compensator  is  thus  implemented  according  to  Equation  (178)  as 


0  10 


0  0  1 


0  0  1 


52  + 


(y  -  yref>  • 


(188) 


The  next  step  in  the  process  is  step  3,  the  design  of  the  full- 
order  stabilizing  observer.  The  structure  of  this  observer  is  given  in  [8]  as 


p  “  (A  -  KjjCJp  +  Kg(y  -  Du)  +  Bu,  u  »  Kp  +  u° 


(189) 


and  it  is  to  be  applied  to  the  plant  so  as  to  provide  a  stable,  satisfactory 
response,  with  Kc  chosen  so  that  (A-KcC)  is  stable.  The  parameter  u°  is  given 
by  Equation  (162).  For  this  example,  the  matrix  (A-Kf.C)  can  be  expanded  to 
give 


(A-KcC) 


0  lj  /kcl\  (1,  0)  T  -kcl  1 

i  —Jl  m 

1  W2/  i-kc2  0 


(190) 


and  the  characteristic  polynomial  of  A  is  found  as 


\  I  -  A 


-1 

“1+kc2  A 


X2  +  kclX  +  (kc2  -  1)  . 


Using  pole  placement  techniques,  place  the  roots  of  Equation  (191)  at 

Xi  -  -4,  \2  “  ~5  , 

which  would  result  in  a  characteristic  polynomial  of 

X2  +  9X  +  20  -  0  . 


(191) 


(192) 


(193) 


By  comparing  Equations  (193)  and  (191),  the  gain  components  can  be  calculated 
to  be 


kcl  “  9»  kc2  "  21  • 

The  observer  given  by  Equation  (189)  can  be  Implemented  as 


\\ .  r-  a 

ty  L'20  °J  \°2 


klPi  +  k2p2. 


(194) 


(195) 


This  observer  was  combined  with  the  plant  model  (with  u°-0)  and  simulated  on  a 
digital  computer.  Appropriate  values  were  chosen  for  the  gain  components  ki, 
k2  to  give  adequate  performance.  These  gains  were  chosen  as 


-40  ,  k2  -  -5 


(196) 


The  last  step  in  the  design  is  step  4  above,  the  evaluation  of  the 
gains  Kq  and  K  for  the  controller  u  and  the  implementation  of  the  complemen¬ 
tary  controller.  The  gains  Kq  and  K  are  to  be  chosen  to  provide  a  satisfactory 
dynamic  response  from  the  augmented  system  described  by  [8] 


A+BK 


8*(C+DK)  -S*DK  C*  \5 


A-KcC  0  II  x-pl  +  0 


(197) 


u  -  (Kg,  0,  K)  I  x-pl, 


(198) 


y  -  (C+DK,  -DK,  0)  lx- pi  +  Du. 


(199) 


These  gains  are  found,  again,  using  pole  placement  techniques.  By  making  the 
appropriate  substitutions  Into  Equation  (193)  an  augmented  system  for  this 
example  Is  expressed  as: 


*2  \  |k01  -39  kQ2“5  40 


k2  k3 


X1“P1 


xi- Pi 


X2“P2 


x2-p2  •  (200) 


Let  the  matrix  In  Equation  (200)  be  denoted  by  A.  The  characteristic  polyno¬ 
mial  of  A  can  be  found  to  be 


X7  "  (ko2  -13)  X6  -  (k01  +8  ko2  -90)  X*  -  (81^  +11  kQ2  -k3  -347)  X* 

-  (llk0i  -20kQ2  +k2  +9k3  -329) X3  -  (-20k01  +kx  +9k2  +20k3  +780 )\2 


-  (9ki  +20k2) X  -  20ki  -  0 


(201) 


If  it  Is  desired  that  the  roots  of  Equations  (201)  be  placed  at 


Xi  -  X2  -  -4,  X3  -  X4  -  -5,  X5  -  X$  -  X7  -  -7, 


(202) 


the  gains  can  be  solved  for  and  the  results  are 


WWW 


Hi* 


KO  “  (*01.  *02)  *  <-348. ,  -26.) 

K  -  <kl,  ^2.  ^3)  -  (-6860., -6027., -2473.) 


(203) 


(204) 


The  robust  controller,  as  given  by  Equation  (162),  is 


Q 


u°  -  (-348.,  -26.)  jxx\  +  (-6860.  ,-6027.  ,-2473.)  /  5l 


(205) 


The  vector  x  is  the  output  of  the  general  complementary  controller  defined  in 
[8]  in  the  form 


x  «  Ax  +  Bu  (206) 
• 

p  -  (A-KcC)p  +  Kc(y-Du)  +  Bu  (207) 
y  «  Cx  +  Du  (208) 
u  «  KP  +  u°  .  (209) 


A  block  diagram  of  the  composite  system,  consisting  of  the  plant 
(Equations  (22)  and  (23)),  the  servo-compensator  (Equation  (188)),  the  full 
order  stabilizing  observer  (Equations  (195)  and  (196)),  the  complementary 
controller  (Equations  (206),  (207),  (208)  and  (209))  and  the  robust  controller 
(Equation  (205)),  is  shorn  in  Figure  27. 

A  listing  of  the  digital  simulation  for  the  system  of  Figure  27  is 
given  in  Table  4.  A  series  of  runs  were  made,  using  this  digital  simulation, 
to  investigate  the  performance  of  the  plant/controller  system  shown  in  Figure 
27.  Figure  28  gives  the  tracking  error  for  cases  with  yref  -  0,  w  -  0  and 
yref  *  0,  w  ■  ec.  Figures  29  and  30  present  the  plant  response  and  tracking 
error,  respectively,  for  a  case  with  yref  ■  1  +  O.lt,  w  -  0.  Figures  31  and 
32  show  the  plant  response  and  tracking  error,  respectively,  for  yref  »  1  + 
O.lt,  with  and  without  the  external  disturbance.  Figure  33  shows  a  comparison 
of  the  tracking  errors  for  four  cases  with  w  -  ec:  (1)  yref  -1+0. It;  (2) 
Yref  *  2t;  (3)  yref  “  -2-4t;  (4)  yref  “  -10+lOt  and  Figure  34  shows  the  plant 
output  response  and  servo-command  input  for  the  case  with  yref  *  -10+10t. 

To  investigate  the  sensitivity  of  the  performance  of  the  control¬ 
ler  design  of  this  section  to  variations  in  the  system  gains,  a  series  of  runs 
were  made  as  discussed  in  the  previous  two  sections.  Figures  35  and  36  show 
the  tracking  error  for  a  case  with  yref  ■  0*  w  -  0  for  +5%,  +10%,  and  +25%, 
+50%  gain  variations,  respectively.  As  shown  in  these  figures,  this  control¬ 
ler  is  not  very  sensitive  to  the  +10%  gain  variations,  however,  the  -50%  vari¬ 
ation  just  about  doubles  the  system  settling  time.  Figure  37  shows  the  re¬ 
sults  of  +25%  and  +50%  gain  variations  when  yref  -  1  +  O.lt,  w  -  ec. 
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TABLE  4.  Simulation  Listing  for  System  of  Figure  27 


COMMON  X<ll),»X<lt>,KUTTA,DT,MX 
DIMENSION  XDAT121) 

NVAR-21 

URITE(20)NVAR 


THIS  PROCRAH  SIMULATES  AH-  EXP  (TIME)  • 
WITH  A  STATE  OBSERVER  FOR  A  SERVO-PROBLEM  • 
FOR  DAVISON'S  METHOD  1  • 


DO  100  1-1,11 
DX(X)-0. 

X(I)-0. 

100  CONTINUE 
X(9)  >  1. 

TIME-0. 

NX-11 

DT-0.01 

PRINT*, '  ENTER  MCO,HCl,rCO,  TCI  ' 
ACCEPT*,  IICO,  WC1,  YCO,  TCI 
XK11  -  -40. 

XK12  -  -S. 

XKC1  «  9. 


XKC2  «  21. 

XKBR  *  0. 

UCC  -  0. 

UOPT  «  0. 

UHAT  «  0. 

U  -  0. 

U  -  0. 

TC  «  0. 

IPRT-0 

1000  CONTINUE 

IF1TIME.GE.10.)  CO  TO  9999 
IPRT-IPRT+1 
DO  200  KUTTA-1,4 
U  -  UC0*EXP(UC1*TIHE) 

YREF  -  YCO  ♦  YC1*TIME 
DX(1)  -  -YREF  ♦  X<1>  ♦  X(8) 

DX(2)  -  X(l) 

DX(3)  -  X(2) 

UOPT— €027. *X(2)-€860.»X<3)-2473. *X( 1 >*UCC 
U  «  UOPT  ♦  XKBR 
DX(4)  -  U  ♦  U  ♦  X(S) 

DX(9)  -  X<4) 


0X<€)  -  XKC1*X(B)  -  XKCt*X<€)  *  X(7) 

DXI7J  -  (l.-XKC2)*X(6)  +  XKCZ*X(9)  ♦  U 
XKBR  -  XKU*X(C>  ♦  XK12»X(7) 

UHAT  -  UOPT  ♦  XK11*X(10)  ♦  XK12*X(11) 

DX(8)  -  UHAT  ♦  X(9) 

DX(9>  -  X(B) 

UCC  -  -2S.*X(8>  -  348.*X(9) 

DX(10)  -  XKC1*X(9)  ♦  X(ll)  -  XKC1*X( 10) 
DX111)  »  (l.-XKC2)»X(10)  ♦  XKC2*X(9)  ♦  UHAT 
CO  TO  (30,60,30,40),KUTTA 
30  CONTINUE 

TIME-TIME*. 5*DT 
40  CONTINUE 
60  CALL  RUNK 
200  CONTINUE 

TRKERR  -  YREF  -  X(9) 

XDAT(1)-X(1) 

XDAT(2)*X(2)  . 

XDAT(3)-X(3) 

XDAT(4)-X<4) 

XDAT(9)-X(9) 

XDAT(6)»X<€) 
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TABLE  4.  Simulation  Listing  for  System  of  Figure  27  -  Continued 


KBAT<7taX<7l 

X9AT<8)-X<8) 

XDAT(9)-X<9> 

XDAT(10)>X<10) 

X9AT(11)*X(11) 

XDAT(1Z)-TIME 

XDAT(13)*DX(1) 

XDAT(14)*U0PT 

XDAT<1S)-U 

XDAT(16)-W 

XDAT(17)*XICBR 

XDAT<18>-TR£F 

XDAT(19)*UHAT 

XDAT(20>-UCC 

XDAT<21>-TRKERR 

IF<IPRT.NE.10)  CO  TO  500 

URITE(20)  <XDAT<I),I-1,NVAR) 

IPRT-0 

PRINT*,'  TINE  «  ',TINE, '  r<T>  -  \X(B>,'  TREF  * 
500  CO  TO  1000 
9999  CONTINUE 
STOP 


EM) 

SUBROUTINE  RUNK 

CONNOR  X(11),DX(11),KUTTA,DT,NX 
DIMENSION  XA<11),DXA(U) 

CO  TO  (10,30,50.70).l0iTTA 
10  DO  20  I-1,MX 
XA(I)«X(I) 

DXA<I)-DT*DX(I> 

20  X(I}-X(I)+.5*DXAU) 

RETURN 

30  TDT-2.»DT 
HDT«.5«0T 
DO  40  I-l.NX 
DXA(I)*DXA(I)+TDT*DX(I) 

40  X(I)-XA(I)«HDT«DX(I) 

RETURN 

90  DO  SO  1*1, NX 
VDT*DT*DX(I> 
DXA(I)*DXA(I)+2.»VDT 
60  X(I)-XA(I)+VDT 
RETURN 

70  DO  80  1-1, NX 


MfREF 


80  X<I)*XA<I)*<DXA<I)*DT»DX<I)>/6 

RETURN 


NG  ERROR 


TRACI 


TIME  <  SEC  > 


Figure  36.  Servo- tracking  error  for  a  case  with  yref  »  0,  vi 
for  ±25%  and  +50%  system  gain  variations. 
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Figure  37.  Servo-tracking  error  for  a  case  with  yref  -1+0 
w  ■  ®c i  y(0)  ■  I  for  +25%  and  +50%  system  gain 


Two  additional  runs  were  made  to  Investigate  the  robustness  of 
this  controller  to  differences  between  the  actual  and  assumed  disturbance  in¬ 
puts  to  the  plant.  Figures  38  and  39,  respectively,  show  the  results  for 
cases  with  w  •  e  and  w  -  e^*'t,  compared  to  a  reference  curve  for  a  case 
with  w  -  ec.  The  results  indicate  that  the  controller  of  this  section  is  less 
sensitive  than  the  controllers  of  the  previous  two  sections  to  these  differ¬ 
ences  . 

The  plots  presented  in  this  section  show  that  this  controller  de¬ 
sign  exhibits  much  larger  initial  transients  than  the  controllers  in  Sections 
V  and  VI.  Apparantly  this  is  due  to  the  influence  of  the  K  term  in  the  equa¬ 
tion  for  u°  (Equations  (162)  and  (205)).  If  this  term  is  removed  from  u°  and 
the  system  is  subjected  to  the  initial  condition  y(0)  ■  1,  the  transient  re¬ 
sponse  is  more  like  that  exhibited  by  the  other  two  controller  designs. 

Figure  40  shows  a  comparison  of  the  system  response,  for  the  design  of  this 
section,  to  the  initial  condition  for  a  case  with  u°  of  Equation  (205)  in¬ 
cluded  and  a  case  with  u°  -  Kqx. 
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Figure  38.  Servo-tracking  error,  yref  *  0,  y(0)  *  1,  for  cases  with: 
(1)  w  -  ec;  (2)  w  ■  e°-5t  . 
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Figure  39.  Servo- tracking  error,  yref  ■  0,  y(0)  *  1,  for  cases 


(2)  w  ■  e^-^t. 
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Figure  40.  Servo-tracking  error  for  a  case  with  yref  ■  0,  w  ■  0, 

y(°)  *  1  and  with:  (1)  u°  as  given  by  Equation  (205) 
(2)u0  with  only  the  Kqx  term  included. 


VIII.  SUMMARY  AND  CONCLUSIONS 

The  controller  design  procedures  given  In  References  [1,2, 3, 8  and  9J 
are  all  straightforward.  In  that,  "design  recipes”  are  given  for  calculating 
the  required  controller  structures.  However,  the  techniques  all  reflect  the 
fact  that  as  the  order  of  the  composite  system  Increases,  the  matrix  manipu¬ 
lations  involved  In  the  design  processes  become  Involved  and  tedious. 

Johnson's  linear  algebraic  design  technique  resulted  in  the  simplest 
controller  structure.  Both  of  Johnson's  design  approaches,  for  the  given  ex¬ 
ample,  resulted  in  an  implementation  with  fewer  additional  integrators  than 
Davison's  technique.  However,  had  the  plant  been  stable,  thus  not  requiring 
the  Inclusion  of  the  stabilizing  observer  in  the  overall  system,  Johnson's 
optimal  controller  and  Davison's  controller  would  have  been  of  about  equal 
complexity.  Johnson's  optimal  controller,  with  the  reduced  order  observer  and 
the  required  optimal  gains,  was  computationally  the  most  complex  of  the  meth¬ 
ods  to  implement. 

As  shown  in  the  plots  in  Sections  V,  VI,  and  VII  and  the  comparison 
plots  in  Appendix  A,  both  of  Johnson's  designs  gave  comparable  performance, 
including  transient  response.  However,  the  design  which  used  the  linear 
algebraic  approach  gave  the  best  overall  performance.  Davison's  design  exhib¬ 
ited  the  worst  transient  behavior,  however,  as  indicated  in  Section  VII,  modi¬ 
fication  of  the  servo-compensator  might  improve  this. 

Johnson's  optimal  controller  design  proved  to  be  the  most  sensitive  to 
variations  in  the  system  gains,  even  to  as  little  as  5  percent  variation.  The 
other  two  controllers  were  able  to  withstand  up  to  50  percent  gain  variations. 
All  of  the  controllers  were  sensitive  to  differences  between  the  actual  and 
modeled  external  disturbance  input. 

When  the  external  disturbance  input  to  the  plant  was  as  modeled,  all  of 
the  controllers  were  able  to  accommodate  the  effects  of  the  external  disturb¬ 
ance  and  achieve  good  servo-tracking  out  to  t  ■  10  seconds.  After  about  10 
seconds,  none  of  the  controller  designs  were  able  to  completely  negate  the 
external  disturbance  effects.  With  no  external  disturbance  included,  all 
three  controllers  exhibited  good  servo-tracking  when  the  simulation  was  run  to 
t  ■  25  seconds,  although  in  some  cases,  the  controller  in  Section  V  began  to 
drift. 
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APPENDIX  A 


To  help  compare  the  performance  of  the  controller  designs  of  Sections  V, 

VI  and  VII,  the  data  shown  on  earlier  plots  in  each  of  those  sections  were 
consolidated.  Each  plot  in  this  Appendix  contains  three  curves  for  each  con¬ 
dition  considered,  one  curve  corresponding  to  each  of  the  controller  designs. 
Figures  A-l  through  A-7  compare  the  servo- tracking  performance  of  the  control¬ 
lers  for  the  various  cases.  Figures  A-8  through  A-13  compare  the  performance 
of  the  controllers  when  the  system  is  subjected  to  gain  variations.  Figures 
A-14  and  A-15  compare  the  performance  of  the  controllers  when  the  system  is 
subjected  to  an  external  disturbance  different  from  the  assumed  disturbance. 

In  addition,  some  results  were  generated  to  demonstrate  the  tracking  per¬ 
formance  of  the  controllers  over  a  longer  period  of  time  than  was  shown  in  the 
previous  plots.  Figure  A-16  is  a  comparison  plot  of  the  system  performance, 
out  to  a  time  of  13  seconds,  when  the  plant  is  subjected  to  an  initial  cond¬ 
ition  of  y(0)  -  1  and  w  ■  ec.  The  performance  of  all  the  controllers,  in  the 
face  of  the  external  disturbance,  is  shown  to  be  good  out  to  about  t  *  12  sec¬ 
onds.  Considering  the  magnitude  of  v  at  t  ■  15  seconds  (w  -  3269020)  the  per¬ 
formance  indicated  is  very  good  but  the  trend  in  the  magnitude  of  the  servo¬ 
tracking  error  is  obvious.  Figures  A-17  through  A-22  show  the  performance  of 
the  controllers,  with  w  -  0,  over  an  Interval  of  25  seconds.  The  designs  of 
Sections  VI  and  VII  performed  well.  The  design  of  Section  V,  as  shown  by  the 
data  in  Figures  A-18,  A-20  and  A-22,  allowed  some  error  buildup. 
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